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With microwave irradiation, the switching current of a Josephson junction coupled to a micro- 
scopic two-level system jumps randomly between two discrete states. We modeled the switching 
process of the coupled system with quantum jump approach that was generally used in quantum 
optics. The parameters that affect the character of the quantum jumps between macroscopic quan- 
tum states are discussed. The results obtained from our theoretical analysis agree well with those of 
the experiments and provide a clear physical picture for the macroscopic quantum jumps in Joseph- 
son junctions coupled with two-level systems. In addition, quantum jumps may serve as a useful 
tool to investigate the microscopic two-level structures in solid-state systems. 
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I. INTRODUCTION 



Recent progress on superconducting qubits based 
upon Josephson Junction (JJ) unambiguously demon- 
strated the quantum behavior of the macroscopic 
variables r^2£££&£ Moreover, quantum jumps, an inter- 
esting quantum phenomenon previously studied in quan- 
tum optics was experimentally demonstrated for 
the first time in a junction coupled with a microscopic 
two-level system (TLS) recently*^ The JJ-TLS coupling 
system possesses A-type energy level structure and mi- 
crowave photons are used to generate transitions between 
quantum states. However, the state of the system is read 
out by detecting macroscopic quantum tunneling process 
rather than that by detecting photon emissions in quan- 
tum optics. Quantum jumps then manifests itself in the 
form of jumping randomly between upper branch and 
lower branch of the switching currents. In the language 
of quantum measurement theory, the switching currents 
in the upper branch or lower branch serve as a pointer 
from which the macroscopic quantum state of the JJ- 
TLS coupling system can be determined. In this situa- 
tion the ensemble description of the dynamics of junc- 
tions based on the master equation metho d 12 ' 13 i 14 fails 
in describing trajectories of a single quantum system. 
Since quantum jump approach developed in the 1980s has 
made great successes in describing fluorescence of single 
trapped ionsi^, in this paper we generalize the quantum 
jump approach to the switching process of JJ-TLS cou- 
pling system and make a systematic study of the param- 
eters that have effects on the process. The same method 
has also been used to investigate quantum jumps in Rabi 
oscillations of a JJ-TLS coupling system^ However, in 
that work the biased current of the junction is fixed at an 
appropriate value while here it keeps changing during the 
switching current measurement. Therefore, new mecha- 
nisms such as Landau-Zener transitions may involve in 
the dynamics of the JJ-TLS coupling system. 



This article is organized as follows. In Sec. II, we de- 
scribe the physics of the current-biased Josephson junc- 
tion briefly and introduce the quantum jump approach 
for simulating the switching process of a current-biased 
junction. In Sec. Ill we generalize the quantum jump 
approach to the switching process of JJ-TLS coupling 
system and discuss the parameters that have effects on 
the process. In Sec. IV we compare our theoretical results 
with experimental data and make a conclusion in Sec.V. 



II. QUANTUM JUMP APPROACH TO 
SWITCHING PROCESS OF A 
CURRENT-BIASED JOSEPHSON JUNCTION 

The Hamiltonian of a current-biased Josephson junc- 
tion as shown in Fig. 1(a) read o 16 ' 17 
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where Iq is the critical current of the Josephson junc- 
tion, / is the bias current, C is the junction capacitance, 
$o = h/2e is the flux quantum, Q denotes the charge 
operator and 5 represents the gauge invariant phase dif- 
ference across the junction, which obeys the convectional 
quantum commutation relation [S, Q] — 2ei. The states 
of the current-biased Josephson junction can be con- 
trolled through the bias current I(t) given by 



I(t) = he + A/(t) = I dc + Ipw COS U)t, 



(2) 



where the classical bias current is parameterized by a dc 
component 1& C and an ac component with the magnitude 
I^u, and frequency to. For Id c < loi the effective po- 
tential of the system (shown in Fig. 1(b)) has a series of 
metastable wells. At low temperature, the current-biased 
junction has quantized energy levels, with the two lowest 
energy states labeled as |0) and |1). Microwaves induce 
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FIG. 1: (a) RCSJ equivalent circuit of a current-biased 
Josephson tunnel junction, (b) Washboard potential of a 
current-biased Josephson junction showing various coherent 
and incoherent processes at low temperature. |0) and |1) are 
ground state and the first excited state which are proposed to 
do quantum information process as a superconducting phase 
qubit. 
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transitions between |0) and |1) at a frequency 
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where w p (J dc ) = 2 1 /4( 2 ^ /$ (7) 1 /2(i _ / dc // ) 1 /4 1S 
the small oscillation frequency at the bottom of the 
washboard potential and AU(Idc) = (2\/2/o$o/37r)(l — 
Idc/Io) 3 ^ 2 is the barrier height. It is apparent from Eq.© 
that the energy spacing u>io is a function of the bias cur- 
rent Idc- Therefore, if we ramp Id c from to Iq, the 
barrier AU is decreasing. At certain Id c called switching 
current the system will tunnel out of the potential well. 
In addition, a microwave with frequency matching the en- 
ergy level spacing will generate a transition between |0) 
and |1). As shown in the top panel of Fig.2(b), the main 
peak of switching current distribution corresponds to the 
tunneling from the ground state |0), and the resonant 
peak corresponds to the tunneling from the first excited 
state |1). By plotting the frequency of microwave vs. the 
position of the resonant peak we can obtain the energy 
spectrum of the junction. 

To simulate the switching process of current-biased 
junction, we firstly write the Hamiltonian of the junc- 
tion in subspace {|0), |1)} 



H 



J J 



Q m cos tot 

fl m COS Lot OJio(Idc) 



(4) 



where Sl m — I^ w y/l/2huJioC is Rabi frequency. Con- 
sidering the dissipative effect of environment, the time 
evolution of the system can be described by the non- 
Hermitian effective Hamiltonian 



ih 

H eff = H qb - y(7io 



ri )|l)<l|-^r |0)<0|, (5) 



where 710 is the energy relaxation rate from |1) to |0), 
and Ti is the tunneling rate from state |i) (i = 0, 1) out of 
the potential (Fig. 1(b)). It is noticed that both 710 and 



FIG. 2: (color online) Simulated switching currents (lower 
panel) of a junction obtained by quantum jump approach 
(a) without microwave and, (b) with microwave respectively. 
The parameters used for simulations are: In = 35.9/iA, 
C = 4pF, oV27r=9.02GHz, and 710=0. 6^s _1 . The ramping 
rate is dld c /dt — 4.5 x 10~ 3 /iA/s. By making the histogram 
of the switching currents we obtained the switching current 
distribution, shown as the symbols in the top panel. The red 
lines are ensemble results obtained using the master equa- 
tions. 



Ti are functions of the bias current Idc 
T, the relaxation rate 710 is given byiS 

^10 Rq 
"2V 



At temperature 



7io = 



R 

where Rq = h/Ae 2 ~ 6.45/cfi is the natural quantum unit 
of resistance and R is the shunting resistance in RCSJ 
model (Fig. 1(a)). The tunneling rate Ti from the state 
\i) can be obtained with the WKB method 
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where T{Ei) is the classical period of motion and Sf(Ei) 
is the action across the classically forbidden region. 

Then the quantum jump approach for simulating the 
switching process of junctions can be summarized as fol- 
lows: 

(i) At t = 0, initializing the junction in the ground 
state: \V(t = 0)) = |0). 

(ii) For I dc (t + At) = I dc (t) + (dI/dt)At, calculate the 
corresponding energy spacing and various transition rates 
according to Eq.([3][7]). 

(iii) Determine whether the system evolves according 
to the schrodinger equation, or makes a 'jump'.— 

(a) If a quantum tunneling escape happens, register 
the switching current / s = Idc(t), and then turn to step 
(v). 
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FIG. 3: (a) (Color online) Schematic of a TLS locating inside 
the Josephson tunnel barrier. Some particles can tunnel be- 
tween two lattice positions with different wavefunctions \L) 
and \R), respectively, (b) Illustration of transition from \g) 
to |e) for the TLS. Assuming the system is initially in state 
\0g), at a certain biased current Id c , the microwave is resonant 
with wio and a transition from \0g) to \ lg) happens. Further- 
more, when the biased current is ramped through the avoided 
energy level crossing, the Landau-Zener transition may lead 
to a finite occupation probability in state |0e). 



(b) If a relaxation event happens, then the system 
jumps to the ground state |0). 

(c) If no jumps happen, the system evolves under the 
influence of the non-Hcrmitian form^ 

(iv) For the case (b) and (c), repeat from step (ii). 

(v) Repeat to obtain the switching current I s . 

(vi) Average switching current I s over many simulation 
runs. 

The numerical results obtained with quantum jump 
approach are shown in Fig. 2. The parameters we used 
in the simulation are from experiments^ In addition, 
we calculate the switching current distribution with the 
master equation method. The agreement between the 
quantum jump approach and the master equation indi- 
cates that quantum jump approach is valid to model the 
switching process of junction. Furthermore, as discussed 
in Sec. Ill, quantum jump approach is more powerful 
than master equation method when stochastic charac- 
teristics of a single quantum system play an important 
role. 



III. QUANTUM JUMP APPROACH TO 
SWITCHING PROCESS OF A JJ-TLS 
COUPLING SYSTEM 

Firstly we give a brief description of the physics of JJ- 
TLS coupling system. TLSs are extensively observed in 
superconducting phas o 18 ' 19 i 20 , charged and Rm£l qubits 
recently. A TLS is understood to be a particle or a small 
group of particles that tunnels between two lattice con- 
figurations, with different wave functions \L) and \R) cor- 
responding to different junction critical current Iql and 
Iqr, respectively (Fig. 3(a)). The interaction Hamilto- 



nian between the junction and TLS can be written as>i£ 

H m = -¥^cob6®\R)(R\-^Zcob6®\L)(L\. (8) 

For convenience, we transfer to the energy eigenstate ba- 
sis of TLS with \g) and |e) being the ground state and the 
excited state, respectively. Then the total Hamiltonian 
of the JJ-TLS in the basis {\0g), \lg), |0e), |le» is given 

by.15,23 

/ Q m cos tot 

H-'h QmCOSUjt UJ W (I dc ) fl c 

fl c ^TLS Q m COS Lot 

V Cl m COSC0t OJio(Idc) + ^TLS 

(9) 

where lotls is the energy frequency of the TLS, and fl c 
is the coupling strength between Josephson junction and 
TLS. In experiments, the coupling strength O c can be 
characterized in spectroscopic measurements and usually 
lies from 20MHz to 200MHz ,II i i8 1 i9 1 20 1 2i 1 22 The time ev0 _ 

lution of the JJ-TLS coupling system under the dissipa- 
tive effect of environments can be described by the effec- 
tive Hamiltonian 

Heff = H- l \ Og \0g)(0g\-^( llo + T lg )\lg)(lg\ 

-yr Oe |0e)(0e| - y( 7l0 +r le )|le)(le|, (10) 

where Ti is the tunneling rate from state \i). We empha- 
size that in the asymmetric double well model of TLS, 
the energy basis of TLS is approximated to the position 
basis. In this approximation, states \g) and |e) corre- 
spond to different critical currents. Therefore, the tun- 
neling rates from different states are different. With no 
loss of generality, suppose the state |e) corresponds to the 
smaller critical current. Then the procedure for simulat- 
ing the switching process of the JJ-TLS coupling system 
can be summarized as follows: 

(i) Initializing the system in state |0<?) for flag = 0, or 
in state |0e) for flag = 1, where flag is a marker. 

(ii) For I dc (t + At) = I dc (t) + (dI/dt)At, calculate 
the corresponding energy spacing uj 10 (I dc) and various 
transition rates. 

(iii) Determine whether the system evolves according 
to the schrodinger equation, or makes a 'jump'. 

(a) If a quantum tunneling event happens, register the 
switching current / s = Idcify- Furthermore, if the system 
tunnels from \0g) or \\g), set flag = 0; else if the system 
tunnels from |0e) or |le), set flag = 1; and then turn to 
step (v). 

(b) If a relaxation event happens, then the system 
jumps to the corresponding ground state |0<?) or |0e), 
i.e., the system jumps from \lg) to |0<ji), or from |le) to 
|0e) . 

(c) If no jumps happen, the system evolves under the 
influence of the non-Hermitian form. 

(iv) For the case (b) and (c), repeat from step (ii). 

(v) Repeat to obtain the switching current I s . 
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FIG. 4: Simulated trajectories of switching current of JJ- 
TLS coupling system under different microwave amplitude. 
The parameters used in the simulations are: u>/27r=9.02GHz, 
lutls/Ztt = 8.7GHz, Q c /2ir = 200MHz, ■y 10 ^0.6/is- 1 , 
dldo/dt = 4.5 x lO^A/s and (a) fi m =2MHz, (b) n m =10MHz. 
With the microwave amplitude increasing , the jumps between 
upper branch and lower branch become more frequent, result- 
ing in a shorter lifetime in each branch. 



FIG. 5: Simulated trajectories of switching current of JJ-TLS 
coupling system under different bias current ramping rates. 
The parameters used in the simulations are: u;/27r=9.02GHz, 
7io=0. 6/is~\ fi m =10MHz and (a) dl dc /dt = 4.5 x 10 3 /iA/s, 
(b) dldc/dt — 8.0 x 10 3 //A/s. With the ramping rate in- 
creasing , the jumps between upper branch and lower branch 
become less frequent, resulting in a longer lifetime in each 
branch. In addition, the switching current in both branches 
becomes higher for the larger ramping rate. 



The simulation results are shown in Fig. 4. It is ap- 
parent that the switching current jumps between upper 
branch and lower branch randomly, which is the major 
characteristic of macroscopic quantum jumps observed 
in experiments^ In addition, it is found that the jumps 
become more frequent with the microwave amplitude in- 
creasing (Fig. 4). The underlying physics can be under- 
stood as follows. The jumps between upper branch and 
lower branch of the switching current are fulfilled through 
the coupling between state \lg) and |0e). With increasing 
the microwave amplitude, the system initialized in |0g) 
has a larger transition rate to \lg) in the expression 



ft 



2(A 2 + 7 2 )- 



(11) 



where 7 = (7 



10 



0g 



Ti g )/2 and A = wio — lj. There- 
fore, it is much easier for the system to jump to state |0e), 
i.e., jump from the upper branch to the lower branch, and 
vice versa. It is easier to understand for the extreme case 
Q m = 0. Then the system has no probability to occupy 
state \lg), thus no probability to transfer to |0e). 

Furthermore, it is noticed that the transition process 
from |1<?) to |0e) is actually a Landau-Zener transition 
as illustrated in Fig. 3(b). Disregarding all decay terms, 
the asymptotic probability of a Landau-Zener transition 



is given b y ' 



LZ 



exp(-27r — £■), 
v 



(12) 



where 2Ml c is the magnitude of the energy splitting, 
and v = de/dt denotes the variation rate of the en- 
ergy spacing for noninteracting levels. Notice that v = 
{de / dldc){dldc/ dt), where de/dldc is determined by the 
intrinsic parameters of the junction and dldc/dt is deter- 
mined by the ramping rate. It can be easily inferred from 
Eq. (|12|) that as the ramping rate increased, the transition 
rate between \lg) and |0e) becomes smaller. Therefore, 
the jumps between upper branch and lower branch of the 
switching currents become less frequent and the lifetime 
for each branch is longer. To support this argument, 
we simulate the trajectories of the switching currents for 
different ramping rates. As shown in Fig. 5, with the 
ramping rate increasing , the jumps between different 
branches become less frequent, as expected from our the- 
oretical analysis. 



IV. EXPERIMENTS 

We have compared the results of our theoretical anal- 
ysis with the experimental data. The sample used in our 



5 




35.70 
35.65 
35.60 
35.55 
35.50 
35.45 



2000 





?•§ '. > 




• : % <|- 






m ;•; 




'; ■ ■'if * 








: •: . I 


.'» ••- *• 

~ 






i = V 


l ' ; 


•j i \ 








t * 







1000 2000 

Number of Switching Events 



< 

3. 



CD 
i_ 

Z! 

o 

o 

- I— ' 

CO 




2000 



35.75 
35.70 
35.65 
35.60 
35.55 
35.50 



MS 



s;. ... .;. '.'' : iv ">'.". - x '/":' ■ 



(b) 



500 1000 1500 2000 

Number of Switching Events 



FIG. 6: Experimental trajectories of switching current of JJ- 
TLS coupling system under different microwave power with 
(a) -7.2dBm and (b) -5.8dBm. The ramping rate is dldc/dt = 
4.5 x 10 3 /iA/s. As increasing the microwave power, the jumps 
become much more frequent. 



FIG. 7: Experimental trajectories of switching current of JJ- 
TLS coupling system under different ramping frequency. The 
microwave power is -6.6dBm and the ramping rates are (a) 
dl dc /dt = 4.5 x 10 3 /zA/s and (b) dl dc /dt = 8.0 x 10 3 /iA/s, 
respectively. With the ramping rate increasing, the jumps 
between different branches become less frequent. 



experiments was a 10/.tmx 10/.tm Nb/AlO^/Nb Joseph- 
son junction. The junction parameters are Iq ~ 36/xA 
and C w 4pF, respectively. The device was thermally 
anchored to the mixing chamber of a dilution refrigera- 
tor with a base temperature of about 18mK. Additional 
three-layer mu-metal surrounding the dewar was used to 
shield the magnetic field. All electrical leads that con- 
nect the junction to room temperature electronics were 
carefully filtered by resistor-capacitor (RC) filters and 
copper powder filters. The center conductor of an open- 
ended coaxial cable was placed above the junction for 
application of microwave. This arrangement resulted in 
>110dB attenuation between the end of the coaxial cable 
and the junction. A saw-tooth bias current was applied 
with a repetition rate of 30Hz ~ 300Hzi n i 26 The junction 
voltage was amplified by a differential amplifier and the 
switching current was recorded when a voltage greater 
than the threshold was first detected during every ramp. 

In the spectroscopy measurement of the junction, an 
avoided crossing caused by the coupling between junction 
and TLS was observed at lo/2tt — 8.7 GHz with an energy 
splitting 2Q c /2ir = 400MHziii When a microwave field 
with oj/2ir = 9.02 GHz was applied, the coupling between 
junction and TLS was turned on. In this case the macro- 
scopic quantum jumps between upper branch and lower 
branch of the switching current were observed. To inves- 
tigate the effect of microwave power as discussed in Sec. 
Ill, we fixed the microwave frequency at u> = 9.02 GHz 
and the ramping rate at dldc/dt = 4.5 x 10 3 fxA/s. The 



microwave power was adjusted from -20dBm to -3dBm. 
As shown in Fig. 6, with the microwave power increasing, 
quantum jumps between different branches become much 
frequent. Similarly, to investigate the effect of ramping 
rate, we adjusted the ramping rate dldc/dt from 2.0 x 10 3 
fiA/s to 16.0 x 10 3 /j,A/s while keeping other parameters 
fixed. As expected, the increasing of ramping rate results 
in less frequent jumps between upper branch and lower 
branch of the switching current (Fig. 7). The agreement 
between our simulated results and the experimental data 
confirmed the validity of the quantum jump approach. 

V. CONCLUSION 

We have used the quantum jump approach to simu- 
late the switching process of the JJ-TLS coupling system. 
The mechanism that dominates the quantum jumps phe- 
nomenon was discussed. In addition, we investigated the 
parameters that have effects on the behavior of quan- 
tum jumps. It is found that a higher microwave power 
or a smaller ramping rate can make the quantum jumps 
happen more frequently, which has significance in con- 
trolling the state of TLS. Furthermore, our theoretical 
results agree with the experimental data, indicating the 
validity of our approach. The model and method we used 
here can be easily generalized to other solid-state systems 
such as flux and charge qubits, quantum dots, trapped 
irons and so on. 
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